On a Condition for Intracellular Adaptive Dynamics for Chemotaxis 
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Microorganisms often perform chemotaxis, (i.e., sensing and moving toward a region with a higher 
concentration of an attractive chemical) by changing the rate of tumbhng for random walk. We 
studied several models with internal adaptive dynamics numerically to examine the validity of the 
condition for chemotaxis proposed by Oosawa and Nakaoka 1], which states that the time scale of 
tumbling frequency is smaller than that of adaptation and greater than that of sensing. Suitably 
renormalizing the timescales showed that the condition holds for a variety of environments and for 
both short- and long-term behavior. 

PACS numbers: 87.17.Jj, 82.39.-k, 05.40.-a 
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I. INTRODUCTION 

Chemotaxis is frequently observed in microorganisms. 
Bacteria move in the direction of a higher concentra- 
tion of an attractive chemical. Although chemotactic be- 
havior is conserved throughout evolution, bacteria never 
move toward a goal directly, but instead often tumble by 
changing direction randomly. By changing the tumbling 
frequency, bacteria assemble around their goal, the re- 
gion rich in attractants. Experimental bacterial systems 
show that bacteria regulate the rate of tumbling and their 
speed, but change direction randomly. Bacteria move 
faster but rarely tumble when located in a region with 
a high concentration of attractant, whereas they tumble 
frequently at low speed in a region with a lower concen- 
tration. However, this experimental observation seems 
peculiar if one considers only normal Brownian move- 
ment, because a high tumbling frequency should yield a 
small bacterial diffusion constant in the continuum limit, 
making the attraction to a high-concentration region im- 
possible. 

By assuming an internal state with adaptation dy- 
namics of a certain timescale into the element exhibit- 
ing Brownian motion, Oosawa and Nakaoka showed more 
than 25 years ago that chemotactic behavior is possi- 
ble 0. In other words, the existence of memory in the 
internal state makes chemotaxis possible. In contrast, de 
Gennes noted recently that the elements can move in 
a favorable direction despite not having an internal state 
to retain some memory. However, this is true only for 
short-term behavior of chemotaxis. As Clark and Grant 
reported , the steady distribution of long-term behav- 
ior is not biased toward a favorable region. A temporal 
change of tumbling frequency with some memory is re- 
quired to have a biased steady distribution. By assum- 
ing a class of autocorrelation function for the change of 
tumbling frequency, Clark and Grant showed that the 
steady distribution of bacteria can be biased toward the 



favorable region. In addition, by introducing some con- 
straints and optimizing both the short-term and long- 
term chemotactic behavior, they obtained a beautiful so- 
lution to the optimal correlation function. 

From a biological viewpoint, the temporal change of 
the tumbling frequency is a result of intracellular dy- 
namics. For any given intracellular dynamics, arbitrary 
autocorrelation function is not possible, and we were in- 
terested in studying the way that specific intracellular 
dynamics allow for chemotactic behavior. Such a condi- 
tion was proposed previously by Oosawa and Nakaoka j3] , 
who studied the steady-state distribution of cells with in- 
ternal adaptive dynamics and the conditions for chemo- 
tactic behavior. When the environment changes, bacteria 
must immediately sense the change and gradually adapt 
to the new environment by exploiting the dynamics of the 
internal state. Oosawa and Nakaoka showed that bacte- 
ria cannot move toward a region with a higher chem- 
ical concentration if the tumbling timescale t is faster 
than the sensing timescale Ts to detect the environmen- 
tal change or slower than the timescale Tq for the adap- 
tation. When tumbling faster, bacteria tumble randomly 
without any directional motion, whereas slower tumbling 
causes the information about direction to disappear be- 
fore the tumbling rate changes. The proposed condition, 
which we call the Oosawa condition, states that the tum- 
bling timescale r is greater than the sensing timescale 
(ts) and smaller than the adaptation timescale (tq). The 
optimal solution by Clark and Grant satisfies this Oo- 
sawa condition. 

In this paper, we introduce a model of an explicit in- 
ternal dynamic system with two degrees of freedom that 
responds and adapts to the external environment. These 
internal dynamics correspond to some intracellular reac- 
tion process, which exhibits both quick-sensing and slow- 
adaptive processes where the parameters Ts and are de- 
rived from the parameters characterizing the intracellular 
reaction dynamics. Using this model, we examined the 
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validity of the Oosawa condition for chemotactic behav- 
ior. We first demonstrate this condition for a chemical 
concentration field with a step function in space. Next, 
using a field with a continuous slope, we describe our ob- 
servations of chemotactic behavior in a broader regime 
than originally proposed. We explain this apparent dis- 
crepancy from the Oosawa condition by the renormal- 
ization of bare timescales through the bacterial motion 
within the slope. Using these renormalized timescales, 
we reconfirm the relevance of the Oosawa condition. 



II. MODEL 

We first introduced an internal state of a cell that re- 
sponds and adapts to the external concentration of the 
signal chemical, and controls the tumbling frequency. 
This internal variable is denoted by u, which might be, 
for example, the concentration of some protein in the cell 
that responds to the external chemical and controls the 
tumbling frequency. This internal chemical responds to 
the concentration of the attractive chemical component 
(termed S here) in the field. As the cell moves and the 
concentration of the external chemical increases, the con- 
centration u increases and returns to the original value, 
a process known as adaptation [Jj. The simplest way to 
have such adaptation dynamics is by introducing another 
internal chemical, whose concentration is given by v, so 
that the change of concentrations is governed by 



dv 



(1) 



We assume here that the fixed point solution u* , v* given 
by f{u*,v*; S) = 0, g{u*,v*) = is stable. If / increases 
with S and u* is independent of S, the response to S 
and adaptation are satisfied because u increases with S 
first, and then returns to the original u*. Here, when 
the solution g{u, v) = involves u but not v, the latter 
constraint is satisfied. For example, g{u,v) = (3uv — jv 
satisfies the condition, where /?, 7 are positive constants. 
In this paper, we report our study of such a case, with 



— = f{u, v; S) ~ S ~ f3uv — au, 
^ = 9{u,v) = I3uv - -fv. 



(2) 



corresponding to the reaction process S* ^ w, m + w — > 2?;, 
and degradation of u and v. Another simple example, 
originally introduced by Othmer Q is given by g{u, v) = 
(S — v)/ fjL, which gives the linear dynamics 



du S ~ {u + v) 

dt rj 

dv S — V 
dt u 



(3) 



with rj and fj, as positive constants. We also briefly dis- 
cuss the result of this model. In both models, after S 
increases, u first increases but later returns to the S- 
independent fixed-point concentration given by u* — ^ / [3 
in the model eq.(2). 

We set the parameter values so that the fixed-point 
solution is stable. In the model eq.(2), this con- 

dition is given by a7 < j3S. In this case, following the 
increase (decrease) of 5", u first increases (decreases) from 
u* and then returns to the original value exponentially 
with time and shows a peak in time. The sensing time 
is estimated as the time to reach this peak and is given by 
^ 1/25*. In contrast, the adaptation time is estimated 
by the relaxation time towards the fixed point, given by 
the eigenvalue of the linearized equation of (2) around 
the fixed point, and is '~ I/7. In this way, the internal 
timescales are represented by the reaction parameters. 

In general, the tumbling frequency is controlled by the 
internal state m, which we assume is given by a continu- 
ous function of the concentration u. The tumbling occurs 
randomly in experiments, whereas the rate changes in re- 
sponse to the external signal 5*. Although the speed could 
also change slightly in bacteria, for simplicity we have 
changed only the rate of tumbling. We assume that the 
cell moves with a constant speed, until it changes direc- 
tion, whose probability (i.e., the rate of tumbling) 1/t{u) 
is given by a function of u. For u — u* , we set the rate as 
l/r* and assume that t{u) is an increasing function for u. 
As an example, we take the form = i-5-tanh(A(-u-Mo)) 

^ ' T\U) T* 

and choose the parameters so that t{u) approaches 2r* 
for u > u* , and t*/2 for u < u* , while A is set sufficiently 
large to respond a change in u by eq.(2). This specific 
choice is not essential and the results we discuss are valid 
provided that the rate (1/t(m)) approaches a value suf- 
ficiently smaller than l/r* for a large u, and sufficiently 
larger than l/r* for small u. This change in the tumbling 
frequency is consistent with experimental data 0. 

Note that the normal tumbling rate 1/t* is indepen- 
dent of the intracellular dynamics, so both are indepen- 
dent of Ts and Ta- This allows us to examine the validity 
of the Oosawa condition. Hereafter, we take the elements 
satisfying eq.(2) and insert them into a one-dimensional 
space, where the external concentration is given as S{x), 
to examine whether the elements assemble in a region 
with larger S{x). 



III. RESULTS ON THE CONDITION FOR 
CHEMOTAXIS 

A. Case 1: Step change in the chemical field 

As a first example, we consider the case with a step- 
function external field (i.e., S{x) = 5+ for a; > and 
S{x) = S~ for X < 0) and examine whether the cells 
assemble in the region x > 0. The number fraction of 
cells at X > 0, numerically obtained for the steady state 
is plotted in Fig. H by changing the parameter r*. In 
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FIG. 1; The number fraction of the ceUs located at a; > (i.e., 
in the S-rich field) are plotted as a function of r * . The fraction 
is obtained from the temporal average for the timespans of 
50,000 (1 < r*), 200,000 (0.1 < r* < 1), and 300,000 (r* < 
0.1), after the steady-state distribution is reached over 100 
cells. See Table|T|for the parameters of the model eq.(2) with 
A, B, and C; the data without internal dynamics are also 
plotted (□) for reference. 
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FIG. 2: The time T to reach a sufficiently large x value, plot- 
ted as a function of r* . The parameter values of intracellular 
dynamics are shown in Table 13 The time T was estimated 
when the center of 100 cells reached x > 1000, after position- 
ing 100 cells initially at s < 10. 



the simulation, all cells were initially located at x < 0. 
The figure shows three examples of different internal 
timescales and r^, and the data without internal dy- 
namics are plotted for reference. Table. ^Jshows the time 
scales(rs , Tq , T\j(to be defined later)) for the three cases 
we used in the plot. 

As shown, the fraction of cells for x > peaks as t* 
changes. For all cases, the most effective tumbling time 
Tpeak that glvcs a peak for the fraction lies at Ts < 
Tpeak ^ t"q. All simulations for other parameters show 
that the elements assemble in the region a; > 0, that is, 
the chemotaxis works well for r, < r* < Tq. These 
results confirm the Oosawa condition. In contrast, as 
T* increases for r* > Tq or decreased for r* < r^, the 
fraction approaches 0.5, implying that chemotaxis is not 
possible. 



TABLE I: The parameters and the timescales. 



Case 


a 


7(= pm 




Ta 


fs 


A 


0.35 


0.25 


0.035 


3 


10 ~ 30 


B 


3.5 


2.5 


0.015 


0.3 


1 ~ 4 


C 


35.0 


25.0 


0.0035 


0.03 


0.2 ~ 0.5 



B. Case 2: Chemical gradient with a constant slope 

Next, we consider the case with an external chemi- 
cal concentration at a constant gradient (i.e., S{x) — 
Sq + sx). Again, by initially positioning all cells in the 
region with small x with low attractant concentration, 
we can quantify whether cells move toward the larger x. 



In this case, the cells with internal adaptive dynamics 
move toward the larger x and stay there, regardless of 
their tumbling time r* . 

The speed for climbing up the slope depends on r* 
and the relation between t* and the internal timescales 
Ts and Ta- To check the speed, we examined the time T 
necessary to reach a specific large x value, as shown in 
Fig. 121 As T* becomes smaller, the time T increases with 
1/t* for T* < T* with some critical value r* that depends 
on Ts and Tq. In other words, efficient adaptive motion 
requires a threshold value r* > t*. The next step was 
to examine whether this range of t* satisfies the Oosawa 
condition given by Tg and Tq. 

We noted that even if r* ^ Ta, the cells can move to 
a larger x efficiently, whereas the critical r* needed to 
increase T 1/r* is much larger than Ts- This seemed 
to violate the original form of the Oosawa condition ( Ts 

< T* < Ta ). 

We note that the actual response and adaptation times 
for cells moving in the described concentration field are 
modified from those estimated from the original intra- 
cellular dynamics. As the cell continuously moves and 
senses the change of the external concentration, these 
timescales 'renormalize'. To examine this effect, we stud- 
ied the intracellular dynamics more closely, because the 
values of Ts and Ta may depend on the dynamics of the in- 
ternal system (u, v) and S is changing continuously with 
time according to the cell's motion. In our model, as 
S{x) continued to change, the equilibrium point of (u, v) 
also changed, which invalidated the baseline u* , v* to give 
Ts and Ta obtained by the linear approximation method 
(Fig. 13- 

Instead, by tracking the time series of u and measur- 
ing the time for the peak and relaxation time to u*, we 
estimated the renormalized values fs and fa. First, as a 
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FIG. 3: The orbit of {u, v) as S changes continuously. The 
large points are the equilibrium points {u*,v*) for given S 
values, whereas the dashed lines are stable manifolds at each 
time. The curved lines with the arrow show the locus of {u, v). 
Before {u, v) reaches the equilibrium point {ul,vl) for S value 
of the time 1, the equilibrium point shifts to (w2, 1^2)- Accord- 
ingly, (u, v) hardly reaches an equilibrium point. 





FIG. 4: The time series of u for the step chemical field S{x) 
(blue) and the gradient with a constant slope (red). The 
series uses parameter values for the case B in Tabled] where 
the fixed-point value u* is 0.5. The response to the change of 
S at t — 0.15 is plotted. 



cell moves to a region with a higher signal concentration 
S, the relaxation to the original u* value hardly occurs 
because, as S increases, the increase in u occurs before 
relaxing to u*. Hence, the renormalized Ta is infinite or 
at least » . Accordingly, one part of the Oosawa con- 
dition T* < Ta is always satisfied. 

The sensing time also changes from Ts when cells are 
placed in a continuously changing field. Although u does 
not relax to the original u*, it increases with S and de- 
creases slowly, as shown in Fig. 0] The renormalized 
sensing time Tg can be estimated by the time needed to 
reach the maximal value. The renormalized fs increases 
from Ts as cells move before u reaches the original peak 
for a given concentration of 5' because the fixed point of 
{u*,v*) continuously changes, as depicted in Fig. |21 

Now, fjj, the threshold beyond which the (renormal- 
ized) Oosawa condition fj, < r* is satisfied, is expected 
to give a criterion for chemotaxis. We compared Ts with 
r* , the critical r* value, to show the increase of T ~ 1/t* 
(see Fig. 01 and found that the two values agree with each 
other. We note that r* > r*(~ fg) is the only con- 
dition for chemotaxis, because T maintains a constant 
value over a wide range for t* > t* . Thus, the (renor- 
malized) Oosawa condition Ts ^ t* < fa correctly es- 
timates the condition for chemotaxis. Finally, we note 
that the renormalized values fs and fa are independent 
of the slope of the concentration of S, and that this con- 
dition gives a criterion for chemotaxis for any gradient in 
the concentration. 

The time T, estimated above as the value necessary 
to reach the region with higher attractant concentration, 
characterizes the ability of chemotaxis over the long term, 
as mentioned by Clark and Grant In contrast, the 
condition for efficient chemotaxis in a shorter-term is dif- 
ferent in general. In our case, however, the chemotactic 



behaviors in the short and long term are not independent, 
but are related through intracellular dynamics. We also 
examined the short-term chemotactic behavior by mea- 
suring the number fraction of cells that moved toward 
the higher attractant concentration at T = 100, starting 
from a random distribution of cells. The fraction is plot- 
ted as a function of t* in Fig. [S] for three examples with 
different internal timescales. Here again, if t* > Ts, the 
fraction is much larger than 0.5, whereas for r* fg, 
the fraction is about 0.5, indicating that the cells move 
randomly. Hence, the Oosawa condition for chemotaxis 
is also valid for short-term behavior. Because the inter- 
nal adaptive dynamics satisfy the renormalized Oosawa 
condition, we conclude that chemotaxis works well over 
both the short and long term. 



IV. DISCUSSION 

We have presented a model of chemotaxis that in- 
cludes internal adaptive dynamics and have shown that 
the chemotactic behavior appears when the Oosawa con- 
dition is satisfied, that is, when the timescale of tumbling 
is greater than the signal-response time and smaller than 
the time for the adaptation. 

Our conclusion about the condition for chemotaxis ap- 
plies generally for a system with intracellular adaptive 
dynamics. For example, we have also studied the lin- 
ear model given by eq.(3) By changing the internal 
timescales Tg and Tq, we have examined the validity of 
the Oosawa condition for chemotaxis. Whether chemo- 
taxis works efficiently is determined by the Oosawa con- 
dition for both the chemical concentration field of a step- 
function and a constant slope, and both for short- and 
long-term behavior. 

To perform chemotaxis, the cell's internal dynamics 
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FIG. 5: The number fraction of cells moving toward the 
larger x, (i.e., for large S{x)) at T = 100. The parameters 
for intracellular dynamics are shown in Table |l] in which all 
conditions are the same as those in Fig|21 The fraction 0.5 
indicates the absence of chemotaxis in short-term behavior. 



have to sense and respond to changes in the the external 
chemical concentration. The timescale of sensing should 
be smaller than that of the tumbling frequency to in- 
duce an effective response. Otherwise, random walks are 
repeated before the response occurs. Thus, the condi- 
tion Ts ^ T* is essential to the short-term response. 
In contrast, response of the tumbling frequency without 
adaptation causes the long-term behavior of cells to be 
represented by a random walk, so that chemotaxis is not 
possible. This implies the need for intracellular adaptive 
dynamics (i.e., relaxation to the original value). How- 
ever, to induce an effective response to external changes, 
the tumbling should occur before the relaxation is com- 
pleted. Hence, the condition r* < Ta is required. Our 
numerical results suggest that the Oosawa condition Ts 
^ T* < Ta is both necessary and sufficient to explain 
the chemotactic behavior. 

Clark and Grant recently reported the conditions 
for the internal response function R{t) to show chemo- 
taxis. In general, the conditions for short- and long-term 
chemotaxis differ, but by imposing some type of opti- 
mization to balance the two behaviors, they were able 
to obtain a certain condition for the response function, 
which implies the existence of a form of memory in the 
autocorrelation function. Our work answers the question 
about why chemotaxis requires dynamics. We can esti- 
mate the response function R{t) from our model, which 
shows that the proper response function in their sense is 
obtained when the Oosawa condition is satisfied. Simi- 
larly, we can define the Oosawa condition in the response 
function R{t) of Clark and Grant. The sensing time Tg is 
estimated from the response function as R{ts) = 0, and 
their solution satisfies Ts t* , whereas is much larger 
than T* . 

Although we recognize the importance of their pa- 



per we believe some controversy remains about the 
optimization condition. We have assumed the existence 
of intracellular dynamics and that these apply to adapta- 
tion after considering the natural properties of intracel- 
lular dynamics j^]. By assuming a class of intracellular 
dynamics, we have shown that chemotaxis occurs in both 
short- and long-term behavior without requiring an opti- 
mization condition, provided that the Oosawa condition 
is satisfied. Because short- and long-term chemotactic 
behaviors are satisfied simultaneously, one expects that 
cells can choose such internal adaptation dynamics that 
enable efficient chemotaxis for a variety of external con- 
ditions. Here we have demonstrated chemotaxis in the 
two cases: with a constant slope and with a stepwise 
change. Because most environmental conditions can be 
represented by the combination of these two cases, the 
Oosawa condition gives a criterion for chemotaxis in gen- 
eral. Furthermore, we have confirmed that chemotaxis 
works when spatial-temporal noise is added in the exter- 
nal environment. 

Because the optimal condition in the internal adap- 
tive dynamics depends on the external motility (tum- 
bling time), the validity of the Oosawa condition can 
be checked experimentally. The table ^ shows the 
timescales estimated from the data obtained for the wild- 
type and mutant cells of E. coli shown in [^0] (see Fig. 3 
and 7 .6]), which supports the validity of the condition. 
In addition, the data in Q suggest that both t* and 
Ta change in mutants whereas the relationship is main- 
tained r* < Ta if the mutant does not lose its chemotac- 
tic ability. In general, different organisms have different 
timescales for tumbling and internal timescales(Ts and 
Ta). It will be important to check the validity of the 
condition in different organisms. 

Living organisms have many reactions with various 
timescales, which, when combined, cause adaptive func- 
tions to emerge. The proper use of different timescales is 
important when generating adaptive output behavior in 
response to changes in external conditions. The process 
of chemotaxis that we have discussed is one of the sim- 
plest mechanisms that takes advantage of the timescale 
differences, whereas the condition for different timescales 
can be generalized in a more complex reaction-network 
system, and in intercellular interactions!^. 

The authors wish to thank Koichi Fujimoto, Shuji Ishi- 
hara, Katsuhiko Sato, Masashi Tachikawa, Satoshi Sawai, 
Naoto Kataoka, and Akihiko Nakajima for useful discus- 
sion and insightful comments. 



TABLE II: The internal time scales of E. coli. (seconds) 



Cell 


T* 




Ta 


Wild- type 


1.5 


0.4 


4 


cheRB"^ 


2.3 


0.6 


3 


cheZ 


12 


2.6 


(larger than 15) 
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